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Abstract 

We show that multiple orthogonal polynomials for r measures (//i, . . . , fir) satisfy 
a system of linear recurrence relations only involving nearest neighbor multi-indices 
fizizej, where ej are the standard unit vectors. The recurrence coefficients are not 
arbitrary but satisfy a system of partial difference equations with boundary values 
given by the recurrence coefficients of the orthogonal polynomials with each of 
measures nj. We show how the Christoffel-Darboux formula for multiple orthogonal 
polynomials can be obtained easily using this information. We give explicit examples 
involving multiple Hermite, Charlier, Laguerre, and Jacobi polynomials. 

1 Introduction 

The three-term recurrence relation is an important piece of information when one is 
studying orthogonal polynomials. If /i is a positive measure on the real line for which all 
the moments 

exist, then we can define orthogonal polynomials {Pn, n = 0,l,2,...}by 

j Pn{x)x^ d^i{x) = A; = 0,l,...,n- 1. (1.1) 

if we normalize the polynomials to be monic, i.e., P„(x) = + - ■ ■ , then these polynomials 
always exist and they are unique whenever the support of /x is an infinite set, and this is 
due to the fact that the Hankel matrix of moments 
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is positive definite for all n > 0. The monic orthogonal polynomials are explicitly given 
by 
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They satisfy a three-term recurrence relation 

^ri+l(a;) = (X - hn)Pn{x) 
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with initial values Pq = 1 and P_i = and recurrence coefficients b„ G R and > 0. An 
important result (the spectral theorem for orthogonal polynomials or Favard's theorem) 
is that any choice of recurrence coefficients 6„ G M and a^+i > (n = 0, 1, 2, . . .) gives a 
sequence of monic orthogonal polynomials for some positive measure n on the real line. 
This measure n is the spectral measure of the Jacobi operator 
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or, when the operator is not self-adjoint, a self-adjoint extension of this operator. For the 
general theory of orthogonal polynomials we refer to Chihara [S], Freud [13], Gautschi 
[13], Ismail [12] and Szego [IH], and for the spectral theory we recommend Simon [18]. 

Recently a generalization of orthogonal polynomials has been considered. Let n = 
{ni,n2, . . . ,nr) G N'' be a multi-index of size |fi| = rii + n2 + . . . + fij. and suppose 
/ii,/i2, ■ ■ ■ , Hr are positive measure on the real line. Type I multiple orthogonals are 
given by the vector (A^ i, . . . , A^ ,,), where A.^ ,,- is a polynomial of degree < rij — 1, for 
which 



0, 



= 0,1, 



\n\ - 2, 
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where n = Hi + ^2 + 
the normalization 



j A^j{x)wj{x) djj,{x) 

+ /ij, and Wj is the Radon- Nikodym derivative dfij/dfi. We use 

r 

^ A^j(x)w;j-(x) dfi{x) = 1. (1.4) 
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The equations fll.3l) - fll.4p are a linear system of |n| equations for the unknown coefficients 
of A.fi^i, . . . , Afi^r- This system has a unique solution if the matrix of mixed moments 
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where 

lyC?) = / x'^dfij{x), j = l,...,r, 



is not singular, in which case we call the multi-index n a normal index. The type II 
multiple orthogonal polynomial is the monic polynomial Pn{x) = x'"' + ■ ■ ■ of degree 
\n\ for which 

J Pfiix)x'^ dni{x) = 0, = 0, 1, . . . ,ni - 1, 

: (1.6) 

Pfi{x)x'' dfir{x) = 0, k = 0,1, ... ,nr — 1. 



This gives a linear system of \n\ equations for the \n\ unknown coefficients of Pfi and 
the matrix of this linear system is the transpose Mi of (11.51) . hence exists and is 
unique whenever n is a normal index. See Aptekarev [2], Coussement and Van Assche 
[20] , Nikishin and Sorokin Chapter 4, §3], Ismail [TS', Chapter 23]. The polynomials 
on the stepline, i.e., 

Pkr+j = P{k+l,...,k+l,k,...,k) 

are also known as rf-orthogonal polynomials (with d = r), see Douak and Maroni [T2] . 
Ben Cheikh and Douak [5|. 

In this paper we assume that all multi-indices are normal and we investigate the 
nearest-neighbor recurrence relations 



xPn{x) — Pfi+ei{x) + hfi^Pn{x) + ^ ] CLnjPn-ej jx). 
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where = (0, . . . , 0, 1, 0, . . . , 0) is the j-th standard unit vector with 1 on the j-th entry, 
and {a^^i, . . . , ,.) and (6,j,i, . . . , bfi^r) are the recurrence coefficients. These recurrence 
relations were derived in [151 Thm. 23.1.11] and it was shown that 



X '■' [X] 

, (1.8) 

/ x"^ ^Pn-s^{x)diij{x) 

and 

hyij = / xPH{x)QH+eAx) dij{x), (1.9) 



where = Yl^i=i ^njWj. There are similar recurrence relations for the type I multiple 



3 



orthogonal polynomials: 



r 



Qri('^) Qn— ei('^) ~l~ ^n— ei,lQn('^) ~l~ '^n jQn+ej (•^) ; 



(1.10) 
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Observe that the same recurrence coefficients (a^_i, . . . , a^j^^) are used but that there is a 
shift in the other recurrence coefficients (^n-e^ j)i<j<r- 

In Section 2 we will show how to find these recurrence relations from a Riemann- 
Hilbert problem for multiple orthogonal polynomials, which was given by Van Assche, 
Geronimo and Kuijlaars |2T]. This approach is not new but merely serves to show that 
these nearest neighbor recurrence relations enter naturally. In Section 3 it will be shown 
that not every choice of recurrence coefficients (a.^,i, . . . , a^i^r) and (6^j,i, . . . , &,j,r) is possi- 
ble: our main result is that the recurrence coefficients satisfy a system of partial difference 
equations. The boundary conditions are given by the recurrence coefficients of the orthog- 
onal polynomials with measure /ij whenever n = nej, i.e., when all the components in the 
multi-index n are zero except for the j-th component. As such the recurrence coefficients 
of multiple orthogonal polynomials are the solution of a discrete integrable system on the 
lattice N''. We will first work with the case r = 2 to simplify the notation and then give 
the general result for r > 2. In Section 4 we will show how the nearest neighbor recurrence 
relations and the partial difference equations for the recurrence coefficients can be used to 
give a simple proof of the Christoffel-Darboux formula for multiple orthogonal polynomi- 
als, which was first given by Daems and Kuijlaars [lU]; see [T] for a Christoffel-Darboux 
formula for multiple orthogonal polynomials of mixed type. This Christoffel-Darboux 
formula plays an important role in the analysis of certain random matrices [6] fl6\ and 
non- intersecting Brownian motions [11] . In Section 5 we will give several examples of mul- 
tiple orthogonal polynomials and give their recurrence coefficients explicitly. For some of 
these examples the recurrence coefficients are given here for the first time. 

2 The Riemann-Hilbert problem 

In this section we will take r = 2 to keep the notation simple. The extension to general r 
is straightforward but requires bigger matrices. Suppose that the measures /xi and /i2 are 
absolutely continuous with weights Wi and W2. If Wi and W2 are Holder continuous, then 
we can formulate the following Riemann-Hilbert problem: find a 3 x 3 matrix function Y 
such that 

1. y is analytic on C \ M, 

2. the limits Y±{x) = lime^o+ Y{x ± ie) exist and 
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3. For z ^ oo one has 




Y{z) = [I + 0{l/z)) \ |. (2.2) 

In case wi or W2 are defined on bounded intervals or semi-infinite intervals one additionally 
needs to specify the behavior of Y near the endpoints of the intervals. In [21] it was shown 
that this Riemann-Hilbert problem has a unique solution in terms of type II multiple 
orthogonal polynomials when {n,m), [n — l,m) and {n,m — 1) are normal indices, i.e., 



C{Pn,mWl) C{Pn,mW2) 
Y = \ Ci{n,m)Pn-l,rn CiC {Pn-l,mWl) CiC{Pn-l,mW2) \ (2.3) 
C2(n,m)P„,^_l C2C{Pn,m-lWl) C2C {Pn,m-lW2) , 



where the Cauchy transform is used 

and the constants ci and C2 are given by 

1 -1 



Pn~l,mix)x'^ Wi{x)dx 



TT- I Pn,m-l(a;)x™ ^W2(x) dx. 



C2{n,m) 27ci 

In case /xi and fi2 are discrete measures on the set A = {xi,X2, X3, . . .} one can formulate a 
Riemann-Hilbert problem for a meromorphic matrix function: find a 3 x 3 matrix function 
Y such that: 

1. Y is meromorphic on C with poles of order one at the points in A, 

2. the residue of Y at x^ is given by 

(0 /ii({xfc}) fi2i{xk}y 





3. Y has the asymptotic condition fl2.2p . 
The solution of this problem is 

\c2{n,m)P^^^.,{z) C2Er=i^^I^/^i({^4) C2Er=i^^^T^/^2({x.}); 

where 

^ CO 

r = ^Pn-l,m{Xk)xl'^ Hl{{Xk}) 



^ 00 

V = y^Pn,m-l{Xk)x'^^^fl2{{Xk}). 



We will derive the recurrence relations for type II multiple orthogonal polynomials from 
the Riemann-Hilbert problem. The result is not limited to Holder continuous weights Wi 
and u'2 but holds in general and can be derived from the biorthogonality [151 Thm. 23.1.6] 

0, if m < n, 
Pn{x)Qrh{w) dfi{x) = \ 0, if \n\ < |m| — 2, 

1, if \m\ = \n\ + 1. 

We merely use the Riemann-Hilbert problem to indicate that the nearest-neighbor recur- 
rence relations come out in a natural way. 

Theorem 2.1. Suppose all multi-indices {n,m) G are normal. Then the type II 
multiple orthogonal polynomials satisfy the system of recurrence relations 

-fra,m+l('^) (-^ dn^rn) Pn,m(^'^) (^n,mPn—l,m(^-^) ^7i,m-Pn,m— 1 (•^) ; (2'5) 

with ao,m = and bnfl = for all n,m > 0. 

Proof. First of all we observe that det Y is an analytic function in C \ M which has no 
jump on the real axis, hence det Y is an entire function. Its behavior near infinity is 
detY{z) = 1 + 0(1/2;) hence by Liouville's theorem we find that dety = 1. We can 
therefore consider the matrix 

Ri{n,m) = Yn+i,mY~^, 

where the subscript (n, m) is used for the Riemann-Hilbert problem with the type II 
multiple orthogonal polynomial Pn,m in the entry of the first row and first column of Yn^m- 
Clearly Ri is an analytic function on C\M, and since l^n,™ and Vn+i,™ have the same jump 
matrix on M we see that Ri has no jump on M. Hence Ri is an entire matrix function. If 
we write the asymptotic condition fl2.2p as 

Y„,uz) = (i+^^ + m/z': 



.n+m 

















where y4(n, m) is the 3x3 matrix coefficient of I/2; in the 0(l/z) term of (12. 2p . then after 
some calculus we find 

' z + Ai^i{n + l,m) - Ai^i{n,m) -Ai^2{n,m) -Ai^'i{n,mf 
i?i(n,m)=| A2,i(n + l,m) \+0{l/z) 

v43,i(n + l,m) 1 

where y4ij(n, m) is the entry on row i and column j of v4(n, m). Liouville's theorem then 
implies that Ri is the matrix polynomial 

' z + Ai^i{n + l,m) - Ai^i{n,m) -Ai^2{n,m) -Ai_3(n, m)^ 
Ri{n,m)=\ A2,i(n + l,m) |, (2.6) 

A3,i(n + l,m) 1 
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and we can therefore write 



In a similar way we also have 



Yn,m+1 = R2in,m)Yn^r 



with 



R2{n,m) 



' z + Ai^i{n,m + 1) - Ai^i{n,m) -Ai^2{n-,m) -^1^3(71, m)^ 
A2,i(n,m + 1) 1 

A3,i(n,m + 1) 



(2.7) 
(2.8) 

(2.9) 



If we now work out the (1, l)-entry of (12. 7p . then we find (12.41) with 

Cn,m = Ai^i{n,m) - Ai^i{n + l,m) 

and 

an,m = ci{n,m)Ai^2{n,m), bn,m = C2{n,m)Ai;s{n,m). 
Similarly, if we work out the (1, l)-entry of (12.81) . then we find (12. 5p with 

dn,m = Ai,i{n, m) - Ai^i{n, m + 1). 



□ 



Observe that deti?i(n, m) = 1 which implies that Ai 2{n,m)A2^i{n + l,m) = 1. If we 
work out the (2, l)-entry of (12. 7p then we find 

Ci(n + l,m)P„,^ = A2,i(n + l,m)P„,^, 

so that ci(n + l,m) = A2,i(n + l,m). Combined with the previous result, this gives 

ci(n, m) 
Ci[n + 1, m) 

which corresponds to (II. 8p when r = 2. Similarly the identity det -R2('^, '^) = 1 implies 
that Ai^3{n,m)A3 i{n,m + 1) = 1 and the (3, l)-entry of (12. 8p gives C2{n,m + 1) = 
A3 i(n, m + 1), so that 

C2(n,m) 



Theorem 2.2. Denote the moments by 



v'l^ = I x^Wj{x) dx 



C2(n,m + 1) 



and the moment matrix by 



( 



(1) ^(1) 
(1) ,,(1) 



(1) „(2) 



K-1 ^' 



(2) 



(1) ^(1) 



^(2) (2) 
'^2n+m-2'^n+in-l '^n+m 



'^n+2m~2/ 



(2.10) 



Then the recurrence coefficients can he written as 

det Mn+i,m det M„_i_^ 



and 



Furthermore we have 



dn,m Cfim 



{det Mn,m) 

det Mn,m+i det Mn,m-l 
( det Mn,m) 

det Mn,m det Mn+i ,m+l 



det M„+i,m det M„ 
Proof. The type II multiple orthogonal polynomial can be written as 



det Mn , 



/I) 



,(1) 



(1) 

1 

(1) 



v. 



(1) ^(2) 
(1) ,,(2) 



,(2) 



/2) 
m—1 
,,(2) 



,(1) 



,(2) „(2) 



(2) ^n+m 



where Mn^m is the moment matrix given in fl2.10p . This formula indeed has all the 
orthogonality properties (II .Op . and from it one easily finds 



x'^Pn,m{x)Wi{x) dx 



:-l)-detM„+i,, 



so that 



In a similar way 



and 



det Mn,ra 

det Mn+i,m det M„ 

— l,m 

(detM„,^)' 

det Mn,m+1 



X"^ Pn,m{x)w2{x) dx 



detM, 



det Mn,m+i det M„,m_i 



( det Mn,m) 

Take the (3, l)-entry of (ET]) to find 

C2{n + 1, m)Pn+i,m-i = As^iin + 1, m)Pn,m + C2{n, m)Pn,m-i- 
On the other hand, if we subtract fl2.4p from (12 .Sp then we have 

If we increase m by one in the penultimate formula, then comparison with the previous 
formula gives 



d. 



C2(n, m + 1) det M„^^ det Mn+i,m+i 



n,m '-^n,m 



C2(n + l,m + l) det M„+i^„ det M„^m+i ■ 



(2.11) 
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Similarly, if we take the (2, l)-entry of (12. 8p . then 

Ci{n,m + l)Pn-l,m+l = A2,i(n,m + l)Pn,m + Ci{n,m)Pn-i,m, 
and if we increase n by one and compare with the previous formulas, then 

ci(n + 1, m) det M^^m det Mn+i,m+i 



ci(n + l,m + l) det M„+i,m det M„,m+i ' 



(2.12) 
□ 



Observe that the arrays c and d can be obtained if one knows the array Ai^i and the 
latter contains the second coefficients of the type II multiple orthogonal polynomials 



Pn,m{x) = + Ai,i(n,m)x 



n+m— 1 



+ 



Each Ai^i{n,m) can also be expressed as the ratio of two moment matrices. 



Ai,i(n,m) 



and the moment matrix in the numerator is 



detM„ 
detM„ 



,,(1) 



V, 



(1) 
1 

(1) 



'^n+m-2 '^n+m-l 
(1) ,,(1) 



V K 



V. 



n+m+1 



(2) 


(2) 



,(2) 



(2) 



2n4 
,(1) 



,(2) 
i+r 
,(2) 



V. 



(2) 

n+m+1 



(2) 



(2) 

^n+2m-3 
'^n+2m-l / 



i.e., it is the moment matrix M„ 
,(1) ,,(1) 



from (12.1 op but with the last row replaced by 



,(1) 



,(2) 



,(2) 



V, 



(2) 

n+2m— 1 



3 Partial difference equations 

We will first consider the case r = 2 to simplify the notation. The case for general r is 
handled in the second part of the section. 

The recurrence coefficients an,m-, bn,m, Cn,m, dn,m are not arbitrary sequences. They 
satisfy certain partial difference equations, which are given in the following theorem. 

Theorem 3.1. Suppose that all the indices (n,m) G are normal. The recurrence coef- 
ficients in the recurrence relations (12. 4 p - (12. 5p for type II multiple orthogonal polynomials 
satisfy the following equations 



dr. 



dn,m 






dg^ 1 '^n+l,r7i dfi 








Cn,m dfi^m 




Cji— l,m dn—l^rn 


h 


r - d 




Cji m— 1 dn m—1 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
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Proof. We obtain these equations by observing that can be obtained in two 

different ways from the information in Yn^m- First we can compute y^+i^m by using (12.41) 
and then we can use (12.51) (with n replaced by + 1) to find y^+i^^+i, i-e. 

Yn+i,m+i = R2{n + 1, m)Ri{n, m)Yn,m- 

On the other hand we can also compute F„,m+i by using (12. 5p and then use (12.41) (with 
m replaced by m + 1) to find 

Yn+i,m+i = Ri{n,m + l)R2{n,m)Yn,m- 

Since det Y = 1 this implies that 

R2{n + 1, m)Ri{n, m) = Ri{n, m + l)R2{n, m). (3.5) 

Observe that 

(z-Cn,m -an,m/ci{n,m) -6n,m/c2(n, m)' 
ci(n + 1, m) 

C2(n + 1, m) 1 

/ z-dn^rn -a^,^/ Ci{n,m) -&„,m/c2(n,m)' 
R2{n,m) = Ci(n, m + 1) 1 

\c2(n,m + l) 



hence (13.51) is equivalent to 



'Si^i{z) 5*1,2(2) ^'1,3(2) 
S2,i{z) 
.^3,1 W 





^0 



















^0 








where 



*S'l,l(2^) i^Z Cn^rn){^Z C^n+l,m) (^n+l,m ^n+l,m (^^ 'J^n,m)(2 Cn,m+l) 

~l~ '^71, m+1 ^n,m+l) 

>-'l,2(^) = ~(^ — "n+l,m) 7 T + (2: — C„,m+1 j 7 r + 



Ci(n,m) ' Ci(n,m) Ci(ri, m + 1) 

•Jl.Sl^;) = — (2; — an+l,mj 7 T 7 — ^ + [Z — Cn^m+1 



C2{n,m) C2(n + l,m) ' C2(n,m) 

5*2,1(2;) = (z - c„.m)ci(ra + l,m + 1) + ci(ra + l,m) - (z - (i„,^)ci(n + l,m + 1) 
5*3,1(2;) = (z - c„,m)c2(n + l,m + 1) - (z - (i„,m)c2(n + l,m + 1) - C2(n,m + 1). 

Observe that 5*1,1 is a polynomial of degree one and the remaining quantities 5*i,2, 5*i,3, 5*2,i 
and 5*3,1 are polynomials of degree zero. If we set the coefficient of z in 5*i,i equal to zero, 
then we find (13. ip . If we put the constant term in 5*i,i equal to zero, then we find (13. 2p . 
If we put 5*1,2 equal to zero and use (I2.12p . then we find (13. 3p . and if we put 5*i,3 equal to 
zero and use ( 12. lip , then we find ( 13. 4p . If we put 52, i and 53,1 equal to zero, then we find 
( 12121) and ( l2lTll . ' ' □ 
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If we introduce the difference operators Ai, A2, Vi and V2 for which 

Aif{n, m) = f{n + 1, m) — f{n, m), A2/(n, m) = f{n, m + 1) — f{n, m) 

Vi/(n, m) = f{n, m) - f{n - 1, m), V2/(n, m) = f{n, m) - f{n, m - 1), 
then we find the partial difference equations 

Aid = A2C, 



Ai(a + 6) - A2(a + 6) 
A2 log a = Vi log(c - d) , 



det 



Aid d 
Aoc c 



Ailog6 = V2log(c-c/). 



The result for general r has more partial difference relations. In fact for each pair 
with i ^ j we have partial difference relations as in Theorem 13.11 for the partial difference 
operators Aj, Aj and Vj, Vj with 

AJ{n) = f{n + el) - f{n), VJ{n) = f{n) - f{n - e,). 

We will not use the Riemann-Hilbert problem explicitly but some of the elements of the 
proof use ideas from the Riemann-Hilbert approach, in particular the transfer matrices 
matrices Rk{n) are inspired by the matrices Ri{n, m) and R2{n, m) which we used for the 
case r = 2. 

Theorem 3.2. Suppose all multi-indices n & W are normal. Suppose 1 < 7^ j < r, 
then the recurrence coefficients for the nearest neighbor recurrence relations ( [i. 7| ) satisfy 



bfi+e^ ,i b. 



n,ii 



k=l 



k=l 



'"it+ei.k 



''n+ej ,i 



det ( ^"+"^^-'' ^ 



(3.6) 
(3.7) 

(3.8) 



Proof. The recurrence relations (11. 7p give a remarkable relation between three neighboring 
polynomials: 

This relation and the k-th relation in (Il.Tp can be written as 



-' n+ej,— ei 



Pn{x) 



\Pn+ek-er{x)/ 

where Bkj{n) 
relation as 



/X — bfi,k — Ctn,l —0'n,2 " " " —Cln,k " " " 

1 Bk,i{n) 
1 Bk,2in) 


















/ Priix) \ 
Pfi—ei (2^) 

[x) 



p 



n-e2 



Pfi-et^ {x) 



VI 0000 Bk,r{n)J \Pn-eM/ 

bn-ej,j — bfi-g.^k (note that Bk,k{n) = 0). We will write this matrix 
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Note that is now a vector containing the multiple orthogonal polynomial and its 
neighbors from below, and not a matrix as in the Riemann-Hilbert problem. Now we 
observe that there are two ways to obtain when i ^ j. The first way is to first 

compute F^+e, from and then to compute from Vri+e,: 

Yfi+e,+ej = Rj{n + ei)Ri{n)Yfi. 

The second way is to first compute Y^j^g^ from Y^ and then to compute YftJ^g-j^g. from 
^ ■ 

^ n+Sj • 

y^+e-+e-. = Ri{n + ej)Rj{n)Yr-i. 
The compatibility between these two ways implies that 

Rj{n + ei)Ri{n)Yfi = Ri{n + ej)Rj{n)Y^. 

The polynomials in the vector Y^ are linearly independent whenever n is a normal index. 
Indeed, if (cq, ci, . . . , c^) are such that 

i=i 

then surely cq = follows by comparing the leading coefficients. If we multiply by x"*"^ 
and integrate with respect to /ijt, then 



SO that Cfc = follows, since the integral does not vanish. Indeed, if the integral vanishes, 
then Pfi — aPfi_g^ is (for every a) a monic polynomial of degree \n\ satisfying the orthog- 
onality relations (11. 6p . which contradicts the normality of n. The linear independence of 
the polynomials in implies that 

Rj{n + ei)Ri{n) = Ri{n + ej)Rj{n). (3.9) 

The transfer matrix Rk{n) can be written as 

- brih -di 



Rk(n) 



BJn] 



where a^j and 1 are (column) vectors of length r containing the a^j- (1 < j < r) and r ones 
respectively, and Bk{n) is a diagonal matrix of size r containing the Bkj{n) (1 < j < r) 
on the diagonal. We then have 

Ri{n + ej)Rj{n) 

[x — bfi^g.^i)(x — bftj) — X]fe=i ^n+ej,k ~ bfi^g.^i)a^ — a^j^g^Bj{n^ 

{x - bri,j)l + B,{n + ej)l -1 ■ + Bi{n + ej)Bj{n) 

and 

Rj{n + ei)Ri{n) 



X — bH+ei,j)ix — bnA) — J2k=l ^n+ei,k —{x — bii+ei,])^^ — di^g,Bi{n) 

(x - brxi)l + Bj {n + ei)l - 1 ■ ai + Bj {n + ei)Bi (n) 
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If we compare the (1, l)-entries in fl3.9p then this gives a polynomial of degree 1 which is 
identically zero. The coefficient of x vanishes whenever (13. 6p holds and the constant term 
vanishes whenever (13.71) holds. If we check the remaining entries on the first row of (13. 9p . 
then we have for 1 < k < r 

bn+ei,jO'n,k ~ Bi^k{n)(ln+ei,k = b^+ej,j0.n^k ~ B j ±{n)a^-\-g. ^k- (3.10) 

If we take this for k = i then we find 

^n+ej,i bfl^g^j bfl^g^^i 

which gives (13. 8p after using (13. 6 p to simplify the denominator on the right hand side. By 
taking = j in (13.101) we find (13. 8p again but with i and j interchanged. □ 

The partial difference relations can be written as 

^ibfij ^jbji^ii 
Aj a^,fe - Ai ari,k = det ( ) , 

k=l k=l OfijJ 

Aj log aH,i = Vj log {bnj - bfi^i) , 

4 The ChristofFel-Darboux formula 

Daems and Kuijlaars [10] have given the following Christoffel-Darboux formula for mul- 
tiple orthogonal polynomials 

Theorem 4.1. Suppose (Tii)j=o,i,...,|n| is a path in N'" such that no = 0, = n and for 
every i G {0, 1, . . . , — 1} one has rij+i — Ui = Ck for some k G {1, 2, . . . , r}. Then 

\n\-l r 

{x-y)^ PnA^)Qri,+i{y) = Pn{x)Qn{y) - afijPH-e,{x)QH+e,{y)- (4.1) 

1=0 j=l 

Observe that the right hand side is independent of the path (nj)j=o,i,- -,|n| ^'^^ ^ 
similar formule for multiple orthogonal polynomials of mixed type is given in [Ij. We 
will give a simple proof based on the nearest neighbor recurrence relations (II. 8p and the 
partial difference equations in Theorem 13.21 

Proof. Use the recurrence relation (II. 7p for Hi and Hi + Ck = Hi+i to find 

r 

xPh, (x) = Ph,+, (x) + bri,,kPn, (x) + ^ an^jPn^^e, (x) . (4.2) 

i=i 

In a similar way we can use the recurrence relation (11.101) for Hi^i and Hi^i — ek = Hi, but 
with X replaced by |/, to find 

r 

yQn,+Ay) = Q^y) + bn,,kQn,+Ay) + J^^^'+u'^^.+i+e.d/)- (4-3) 
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Multiply fl4.2p by Qui+i iy) and fl4.3p by P^. (x) and subtract both equations, to find 

r r 

+ ^ari^jPrt^_g^{x)Qa^^^{y) - ^ a^,^, jP^^(x)Q^^^,+e-- (y). (4.4) 
i=i j=i 

If we sum this for i = to |n| — 1, then the left hand side corresponds to the left hand 
side in (14.11) . and the first terms on the right hand side give 

|ra|-l 
1=0 

by using the telescoping property and = 0. The two sums on the right hand side of 
(14. 4 p need to be rewritten a bit. Observe that the recurrence relation (II. 7p implies that 

If we use this for n = Hi — ej, then this gives 

PriXx) = Pn,+ ,-e,{x) + (^n.-e, ,fc - h^.^g^ j)Pfi^.s^{x) . (4.5) 

In a similar way, the recurrence relation (ll.lOp implies 

Qn-eAv) - Qn-e,iy) = (&n-e,j " &ri-efe,fc)Qn(l/) • 

If we use this for n = n^+i + ej, then we find 

<5n.+l (y) = Qn,+e, (Z/) + (&n,+e,,fc " &n,+ij)Qn,+i+e- (Z/)- (4.6) 

Replace P^^ (x) in the second sum of (14. 4 p by (14. 5 p and replace Qn,+i (y) in the first sum 
of (031) by (USD, then we find that 

r r 
^ ^ni,jPni-ej{x)Qni+i{y) ~ ^ (^ni+i,jPni{x)Qni+i+ej jy) 



r r 



^ ^ ^nijPfii-ej ix)Qni+ej (l/) ^ ^ (^Ui+ijPfii+i-ej {x)Qni+i+ej (l/) 
r 

+ Pn—Cj (3;)Qni+i+ej (y) [o^. (6^.+^^. ^fc — 6ni+i,j) ~ O-Hi+ijibfii-ej ,k " ^n— e^-j')] • 

i=i 

The last sum on the right hand side vanishes because of the partial difference relations in 
Theorem 13. 2[ in particular (13. 8p and (13. 6 p but with i = k. The other two sums are now 
of the form Si — S'j+i, so that we are summing a telescoping series, giving 

|n|-l r 

^ ^ (Si — Si^i) = Sq — S*!^! = — ^ ^ O'njPft-ei {x)Qn+ej (Z/) ; 
i=0 3=1 

which is the sum on the right hand side of (14. ip . □ 
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5 Some examples 



In this section we give several examples of multiple orthogonal polynomials and their 
nearest neighbor recurrence coefficients. For each example one can verify that the partial 
difference equations in Theorem 13.21 indeed hold. Some of these examples were given 
before, e.g., in [3], [1], |5], |7], [9], [151 Chap. 23], |20], but the explicit expression for the 
recurrence coefficients appears here for the first time for most of these families. 

5.1 Multiple Hermite polynomials 

See [201 §3.4], ^ §23.5], [7]. These are given by 

oo 

x'=i/^(x)e-^'+'^^-^ dx = 0, k = 0,l,...,nj-l, 

■oo 

for 1 < J < r, where Cj ^ Cj whenever i j. The recurrence relation is explicitly given as 



Cfc 1 

xH^ix) = Hfi+g^{x) + —Hfi{x) + :^ njH, 
for 1 < A; < r, so that 



^n,j = Cj/2, anj = nj/2, 1 < j < r. 

5.2 Multiple Charlier polynomials 

See [H §4.1], [151 §23.6.1]. The orthogonality relations are 

oo 

5^Q(A;)A;^^ = 0, £ = 0,l,...,n,-l. 



fc=0 



for 1 < j < r, where > and Oj ^ aj whenever i ^ j. The recurrence relation is given 
by 

r 

xCft{x) = Cr-i+e^ix) + (cfc + |n|)C^(a;) + '^njajCH-g^{x), 

i=i 

so that 

bnj = \n\ + aj, aaj = rijaj, I < j < r. 

5.3 Multiple Laguerre polynomials of the first kind 

See [20l §3.2] and [151 §23.4.1], [7]. These are given by the orthogonality relations 



oo 



x^LH{x)x°''e dx = 0, A; = 0, 1, . . . , rij - 1, 



for I < j < r, where ai, . . . , > —1 and ctj — aj ^ Z. They can be obtained using the 
Rodrigues formula 

^ ( d'^i \ 
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where the product of the differential operators can be taken in any order. This Rodrigues 
formula is useful for computing the recurrence coefficients. Indeed, we easily find 

roo POO jn\ / Jrii \ 

/ a;"^L^(x)a;"^e-" = (-1)1"! / TT a;""'- — e"" rfx 

7o Wo dx^^ V dx^^ ) 

and integration by parts {rii times) gives 

r°° / (jrii \ _!_ f AUi \ 

^ ' Jo \dx^^ J Ail^ dx". J 



Repeating this r times gives 



/OO ^ / \ \ /*oo 

= r(n,+a, + l)n ' ' W- 
j=i ^ * ^ 



If we then use (11. 8p then we find 

r 

rii + a, — a 



n — — j=i,---,^- (5.2) 



The recurrence coefficients 6^,^ can be obtained by comparing the coefficients of x'"' in 
the recurrence relation (II. 7p . They are 

&n,fc = 1^1 + ?^fc + a/c + 1, A; = l,...,r. (5.3) 

These recurrence coefficients are here given for the first time. Observe that for r = 1 
one retrieves the well known recurrence coefficients for the monic Laguerre polynomials 
= n{n + a) and 6„ = 2n + a + 1. 

An interesting observation is that these recurrence coefficients of multiple Laguerre 
polynomials of the first kind satisfy 



r / r 



i=i i=i 
Indeed, if we introduce the polynomials 



i=i j=i 
then a^j = (ra^ + aj)qr{nj + Oij) /Q'^,^{nj + a^) and hence by the residue theorem 

zqr{z) 



5^"- = 2^/ 

3 = 1 



dz 
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where F is a closed contour around all the zeros + a, {i = 1, . . . ,r) of Qr,n- We can 
take for instance the circle = {i?e*^ : 6 G [0,27r]} for R large enough. If we change 
the variable z = 1/^, then the integral becomes 



where Cr is the circle with radius r = 1/R and the orientation is positive by changing the 
sign appropriately. Here q*{z) = z^qr{l/z) and Q*.^{z) = z'^Qr^ni)-/ z) are the polynomials 
taking in reversed order. The residue theorem now tells us that the integral is the residue 
of IriO/C^QlniO at ^ = 0, and this can easily be computed by using 

1 — fry- 



1 - ^{rij + aj) 
so that 



and the residue of (IriO/^'^QrniO thus corresponds to the coefficient of the quadratic 
term, giving ( 15. 4p . 



5.4 Multiple Laguerre polynomials of the second kind 

See [m Remark 5 on p. 160], [20, §3.3] and [15, §23.4.2], [7]. These are given by the 
orthogonality relations 

oo 

x''Lfi{x)x°e~''^'' dx = 0, k = 0,1, ... ,nj — 1, 

for 1 < j < r, where a > —1, Ci, . . . ,Cr > and Cj ^ cj whenever i j. They can be 
obtained using the Rodrigues formula 

j x'^L^x) = f[ (e'^'^^e-^^"^ (5.5) 

where the differential operators in the product can be taken in any order. A useful integral 
is 

which can be evaluated by using integration by parts in a similar way as in the previous 
example. Observe that the right hand side has a zero at A = Cj of multiplicity rij. Using 
(15. 5p we thus have for A > 



Clearly 

d^ /-oo 



dA^ 



e-^^a;"L„^(x) dx 



oo 

= i-l)^ I x^e'''''x''Ln{x) dx = 0, 0<k< n^, 

A=c, 
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which confirms the orthogonahty relations, and for k = Uj 
If we then use (11 .Sp then we find 

(\n\ + a)ni 

aH,j = ^' ' 2 , 1 < J < ^- (5.6) 

For the coefficients bfi^k we compare the coefficients of x'"' on both sides of the recurrence 
relation f ll.7p and use the exphcit expression (23.4.5) in [TJ] to find 

I — # I ^ 

72 + « + 1 v-^ IT- 



Ck 

J 



+j:x (5.7) 



Observe that for r = 1 and Ci = 1 we retrieve the recurrence coefficients of monic Laguerre 
polynomials bn = 2n + a + 1 and = n{n + a). 

5.5 Jacobi-Piiieiro polynomials 

See [TTl Remark 7 on p. 162], [201 §2-1] and [151 §23.3.2]. These are multiple orthogonal 
polynomials on [0,1] for the Jacobi weights xf{l — x)^, with ai, . . . ,ar, /3 > —1 and 
ctj — aj ^ Z. They satisfy 

[ Pn{x)x^x'^^ {l-xfdx = Q, A; = 0, 1, ... , - 1, 1 < j < r. 
They are given by the Rodrigues formula 

r 

(-1)1"! JJd^l + + /3 + l)n, (1 - x)''P„Kx) 

=n(^""^^^"^^"0(i-^)"""'' (5-8) 

where the product of differential operators is the same as for the multiple Laguerre poly- 
nomials of the first kind. One has 

Jo nLid^l + «i + + Ik r(|n|+/? + 7 + 2) 

(see [151 §23/3.2] which has an extra (— l)'"' that shouldn't be there) so that 

f _ - r*r-)".! r(„, + a, + i)r(|S| + + 1) 



riLid^l + ai + P + 1)„, r(|n| + rij + aj + (3 + 2) 
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If we use fll.81) we then find for 1 < j < r 



TT "-j ^ - TT 

. ^\^.nj - ni + aj - ai ^\ \n\ +ni + ai + 13 







(|n| +nj + + 




+ /3)(|n| 





We can write this as 

qr{-\n\ - /3) + oij) 



X 







1 




(1^1 


+ nj +aj + l)(|n| 


+ nj + a 




+ n^- + + /3 - 1) 



where the polynomials and Qr,^ are given, as before, by 

r r 

Qr{x) = Y\_ix - aj), Qr,n(a:) = - - ftj). (5.10) 

i=i i=i 

Observe that for r = 1 we retrieve the recurrence coefficients 

2 n{n + a){n + f3){n + a + (3) 



{2n + a + 13 -l){2n + a + l3y{2n + a + /3 + l) 



for monic Jacobi polynomials on [0, 1]. 

The other recurrence coefficients can be obtained by 6^^^ = 5.^ — 5n+4 , where 5^ is the 
coefficient of x'"'^^ of the polynomial Ph{x): 

Pn{x) = + + • • • . 

For Jacobi-Pineiro polynomials we have 

where the polynomials and Qr,™ are given in (15.101) . To see that this is true, one can 
use the Rodrigues formula (15. 8p to find that 

(x°'+i(l - xf'-'Pttf-^'-'\x)y = -i\n\ + a, + /3 + l)a;"»(l - xfpjf^^\x) 

(see also Eq. (1.4)]). If we identify the coefficients in this identity, then 

{\n\ +ai + f3 + l)Sf{{d, (3) = (|n| + + (3+)6H-gX^ + e^, /3 + 1) - \n\ - ai. 

This can be used to prove (15. lip by induction on r. Observe that we can use a decompo- 
sition in partial fractions to write 



3=1 
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so that 



and 



qr{-\n\ -/3) ^"^^^ ^ "i)/'5r,nK- + «j) 



Qr,n(-|^| - /9) ^ |n| + + q;j + /3 



Off ^ -\n\ + {\n\ + /3) - 



3=1 



\n\ + rij + aj + P 



Furthermore 



\n\ 



so that a httle calculus gives 

|fi| + rij + + ^ 

Note that for r = 1 this indeed gives 

_ —n{n + a) 

from which 

On — (^n — (^n+l 



2 2(2n + a + /3)(2n + a + ;g + 2)' 

which is indeed the recurrence coefficient 6„ for the Jacobi polynomials on [0, 1]. 
We can find a nice closed form for the sum Y7j=i ^n,j- We have 

n-i ^ R\ ^r{-\n\ - /3) 
2_^an,j = (|n| 



Qr,n{-\^\ - 1^) 



<3r,nK- + "i) (1^1 + % + + /5 + + + aj + /3)(|n| + rij + aj + l3-l)' 

and for the sum in this expression, we can use the residue theorem to find 

^ qr{nj + a,) nj + a,- 



E 



I Q'rfii'^j + «j) (1^1 + + Qij + /3 + l)(|n| + Uj + aj + /3){\n\ + Uj + a^- + /3 - 1) 

1 f qr{z) Z dz 



2m Jr Qr,n{^) {z + \n\+P + l){z + \n\ + P){z + \n\ + P - 1)' 



where F is a closed contour around all the zeros rii + ai {i — 1, . . . , r) of Qr-.n but with 
— \n\ — /3, —\n\ — /3 ± 1 outside the contour. If we choose F so that is inside and change 
the variable z — 1/^, then the integral becomes 



27ri Jr 



r* QUO (1 + m +p+mi+ m + +m+p- 1)) ' 
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where F* is now a contour around with the poles — l/(|n + /3) and — l/(|n| +/3± 1) inside 
and the zeros of Q*^ outside. The polynomials q* and Q*^ are the reversed polynomials 
(as in §5.3p . The residue theorem therefore evaluates this contour integral as 

1 qr{-\n\ - /3 - 1) ^1^1 , ^ g^(-|r?|-/3) ^i^,,^^,! qr{-\n\ - P + I) n^i , p , ,\ 

2Q,«(-|S|-fl-l)(l"l+^-l'-Q.,.,(-|B|-fl)<l"l+«+2Q„,(-|B|-^ + l)'l"|-'''+''' 

and we conclude that 

g = + Q-FN^ U a,.(-N-/3-i) ^'"' + ^ - 

~ -(|n| + /3) + 1 ^-^;'"l~V^\ (|n| + /3 + 1; 



Qr,n{-\n\ - 13) 2Qr,n{-\n\-f3 + l) 

Observe that the expression between brackets is —1/2 times the second difference AV of 
zqr{z)/Qr,n{z) evaluated at z = —\n\ — /3 and that this expression vanishes when \n\ = 0. 
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